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Abstract 

In super symmetric models with R parity violation, constraints on 
superpartner masses are significantly weaker than in models which 
conserve R parity. We find in regions of parameter space where a 
neutral gaugino or third generation scalar is decoupled from the Z 
and/or traditional gaugino mass relationships do not hold allow some 
particles to be light enough to allow them to be the decay products of 
the Higgs. For example, surprisingly a stau could be lighter than 30 
GeV if it decays hadronically. We estimate the Higgs bounds when one 
of these decays dominates and find allowed Higgs masses well below 
the current LEP bound on the standard model Higgs. We also survey 
the rich variety of final states in Higgs production. 
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1 Introduction 



Softly broken supersymmetry (SUSY) is arguably our best candidate for an 
explanation of what cuts off the quadratic divergent quantum corrections to 
the squared Higgs mass in the standard model (SM). Because of these correc- 
tions, the SM alone, with a cutoff much higher than the weak scale, suffers 
from a fine-tuning problem. The mass scale of superpartners effectively rep- 
resent the cutoff and above that scale the tuning is no longer necessary. Thus, 
to completely remove tuning from the Higgs, the natural scale for supcrpart- 
ner masses is on the order of the weak scale. 

Of course neither the Higgs superpartners have been discovered. The 
bound on the SM Higgs mass is 114.4 GeV [?], and this bound, to a good 
approximation, applies to the lightest Higgs in most of the parameter space 
of the so-called 'mSUGRA' version of the minimal supersymmetric standard 
model (MSSM) [?]. These bounds are strong enough to require fine tuning 
within the mSUGRA model, the most restrictive being the Higgs mass bound 

- the physical Higgs mass in mSUGRA lies naturally below the Z mass, and 
increases above the Z with contributions only logarithmically sensitive to the 
superpartner mass scale. Superpartner masses (mainly the stop mass) need 
to be pushed beyond the electroweak scale in order to satisfy the bound. 

Going beyond mSUGRA, Higgs and superpartner mass bounds can change 

- and sometimes radically. Here we study a version of the MSSM with R- 
parity violation. We also consider regions where the standard gaugino mass 
relations do not hold. The main phenomenological feature of violating R- 
parity is that the lightest superpartner is no longer stable and that 'missing 
energy' is either strongly reduced or eliminated in events with superpartners 
produced. This single change, allowing the lightest supersymmetric parti- 
cle (LSP) to decay, reduces bounds on nearly every superpartner (save the 
chargino) to below 100 GeV, and sometimes well below. This keeps open the 
possibility of the Higgs boson decaying into superpartners. 

In this paper, we explore Higgs decays into a pair of LSPs which result in 
multi-body final states. If the Higgs decays in a non-SM way, it affects the 
lower bound on the Higgs mass as the standard LEP searches are less efficient 
or simply do not apply. Based on a wide range of searches, we estimate that 
many of these new decays reduce the lower bound on the Higgs mass to 105 
GeV or less. While this seems like a small change, the exponential sensitivity 
of the superpartner (e.*?., stop) masses to the needed Higgs mass correction 
makes this reduction relevant for re-opening SUSY parameter space. The 
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key is beating the standard decays of /i ^ 66 by a large amount - and 
this happens in broad regions of parameter space the Higgs sector is in the 
decouphng hmit. 

This paper is organized as follows, section 2 reviews R parity violating 
operators and their bounds. Section 3 discusses constraints, signatures and 
parameter space for neutral gaugino LSPs. Section 4 discusses constraints, 
signatures and parameter space for third generation scalar LSPs, and section 
5 concludes. In the appendix, we discuss the fit to e^e~ hadronic cross section 
data when a new light charged scalar is in the spectrum. 

2 R parity violation 

In the MSSM generally one imposes a symmetry known as R parity whereby 
fields and superfields are given opposite parity. Only operators with positive 
parity are allowed to appear. The LSP is stable; if one starts with N super- 
particles, one must end with N mod 2 due to their negative parity. However 
we are free to introduce parity violating operators into the superpotential: 

W D fiiLiH + XijkLiLjEl + X[jf^LiQjDl 

+ K^kUWPl (2.1) 

where L, E'^, H, Q, W^, and D'^ are lepton doublet, lepton singlet, up-type 
Higgs, quark doublet, up-type quark singlet, and down-type quark singlet 
superfields respectively, and the ijk are flavor indices. The first operators 
violate lepton number conservation, while the third violates baryon num- 
ber conservation. An acceptable theory cannot have both LNV and BNV 
couplings of non-negligible size as it would cause rapid proton decay. 

For a general flavor of RP V bounds see [21 S] ■ In general LLE operators 
are of order a few times 10~^, while LQD and UDD operators are bounded at 
10^^ though a few couplings have much harsher bounds and some are order 
1. LNV bounds come from a great variety of sources. Semi-leptonic meson 
decays put bounds on a many products of couplings [5] most of which fall in 
the range 10~^ — 10~^. The most stringent bound for LLE operators is on 
the Ai33 operator and comes from contributions to the electron neutrino' ma- 
jorana mass and is less than 10~^ for 100 GeV superpartner masses. Other 
LLE bounds come from measurements of Rr and Rrn- The LQD constraints 
come from a variety of sources. Neutrinoless double beta decay puts a strict 
bound on the \[^^ operator of ~ 10"^ x {m^/lOOGeVY^^im^/lOOGeVy. 
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Charged current universality, Dg and D meson decays and the measurement 
of Rt- also place bounds and in general operators involving third generation 
couplings A33]^, A332, A333 and A232 are relatively unconstrained. The strictest 
BNV bounds come from double nucleon decay and apply to first genera- 
tion couplings A^^2 ~ 10~^. Neutron anti neutron oscillation limits A^]^3 to 
~ 10~^. Other couplings are less constrained and bounds vary from 10~^ to 
1. A generally safe statement is that with the exception of a few pathological 
couplings cited in the canon, nearly all bounds allow for reasonably prompt 
decay of a superpartner into SM particles. Avoiding stringent constraints on 
some of the first generation RPV operators may be remedied easily by requir- 
ing a Yukawa-hierarchical coupling scheme where third generation couplings 
are dominant; this scheme will be of particular value when considering LNV 
decays of gaugino LSPs. Here we favor a Yukawa hierarchical scheme and 
consider one only one RPV coupling at a time to be turned on. We now go 
on to list possible LSPs which both may be hghter than half the Higgs mass, 
and then may decay via RPV. 

3 Gaugino LSPs 

3.1 Topology and General Constraints 

The Higgs may decay into a pair of neutral gaugino LSPs. The RPV decay 
then proceeds as each gaugino decays to a fermion and virtual sfermion which 
itself decays via an RPV vertex to two SM particles. It would be simple here 
to think of gauginos decaying through an effective four fermi interaction to 
three fermions Xodi^jk/Mj. The topology of these decays is Higgs to 6 SM 
fermions. In the case of BNV we would have Higgs to six jets. For LNV, 
Higgs to 4 jets plus 2 leptons, and LLE to 6 leptons. 

There are several general constraints on this scenario applicable to all 
gaugino LSPs. The first is the mass requirement that gauginos be less than 
half of the Higgs mass. What usually stands in the way of gauginos this 
light are assumptions about mass relations in standard SUSY mass schemes. 
Gauge mediation, anomaly mediation, and MSUGRA relate the gauginos 
by a single mass parameter. There is a tight experimental lower bound 
on the chargino mass of 102.7 GeV, which then constrains the mass of the 
lightest neutralino through standard gaugino mass relations. In addition, 
most standard mass schemes also predict a gluino much heavier than the 
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other gauginos. For example, MSUGRA predicts the M2:M3 ratio to be 
about 3:1 thus a chargino passing the lower mass bound implies a heavy 
gluino. By rejecting the standard gaugino mass relationships we open up 
large regions of SUSY parameter space. 

We must also address lower mass bounds on sparticles from existing RPV 
searches. DELPHI and L3 have put lower bounds on the neutralino as the 
LSP decaying through RPV [S] [Z] [B]- L3's mass limits for xo decaying via 
LQD is 32.6GeV, while it's limits on decay via LIE and UDD are 40.2GeV 
and 39.9GeV respectively. DELPHI has put lower mass limits on decay 
via LLE and UDD as 39.5GeV and 38GeV. At DELPHI, a simple counting 
experiment was done assuming gaugino mass unification. We will examine 
the LLE decay count, since it was the lowest; if we avoid the LLE constraint 
we avoid them all. In the LLE channel DELPHI found only 1.5 events above 
the background at 95 percent confidence. Recalling that = eLa where 
e is the efficiency, we see that for e between .11 and .38 and an integrated 
luminosity of 437.8p6~^, the production limit is ~ .03-. 01 pb. If we assume 
no contribution from charginos, a selectron mass of between 380 and 500 
GeV will sufficiently suppress this process for a 30 GeV neutralino. There 
were additional searches at ALEPH that placed lower bounds of 29 GeV and 
23 GeV on LQD and LLE decays of neutrahnos respectively [9] [1^. These 
analyses also relied on MSUGRA. For example, they put an upper bound 
on neutralino production from LLE decay which is at least .5pb. Without 
MSUGRA, we may easily beat this production bound by picking the selectron 
mass. 

In the case of gluino as the LSP no direct decay bounds are quoted. An 
indirect search exists which looks for direct decay of a neutralino produced 
in the decay of a squark, and assuming gaugino mass unification [TT]. This 
search concerns the coupling Ag^^ only and looks for muons in the final sate. 
For reason which will be explained this does not fall into our scenarios. In- 
terestingly, a lower mass limit on the gluino is set at 6.3 GeV [12]. This 
is obtained from the contribution of hadronic decay width of the Z when 
e+e^ qq and a quark radiates a gluon which decays to a gluino pair. 

The decay of Higgs to neutralinos must also beat the standard Higgs 
decays if current searches are not to have already ruled out the Higgs. The 
largest standard decay rate is Higgs^ bb. The bottom Yukawa coupling is 
not particularly large and we shall show that it does not require special fine 
tuning in parameter space in order for Higgs decays to gauginos to beat the 
bottom decay by a factor of 5 or so. 
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The next constraints heavily involve the couplings to Z bosons. Decays 
of light neutral gauginos may not substantially change the Z width, nor may 
they have a large effect on the predictions for total e+ e- to hadrons cross 
section. These issues may be solved by decoupling the lightest neutralino 
sufficiently from the Z, and making sure loop induced couplings from of the 
Z to the gluino are not too great. 

There is a constraint from supersymmetric contributions to 6 — > 57 which 
come mainly from the charged Higgs and the chargino (diagrams involving 
the gluino also make contributions but are suppressed by insertions of the 
CKM matrix). Diagrams with the ^ind have a stop in the loop. 
Avoiding large contributions to this process is accomplished in two ways: if 
the stop and the charged Higgs are heavy and all diagrams are suppressed, 
or if the diagrams cancel. Often there is some tension between a heavy stop, 
which could tune the Higgs sector, and small b — sj. The current experimen- 
tal bound is 3.55 + -.2Atf + -.03x10-^^. The current Standard Model 
theoretical prediction made at NNLO is now significantly below the measure- 
ment at 3.15H — .23x10-^(11], so there is some small space for supersymmetric 
contributions. 

There are generational constraints on LNV decays which have charged 
leptons in the final state. Tight constraints from Tevatron exist on like sign 
dilepton signals [15]. This search found a slight excess of like sign dilepton 
events, 20 events at 95 percent confidence, at an integrated luminosity of 
lfb~^. The efficiencies for seeing new physics that produced a WZ were as- 
sumed to be around 8 percent. In SUSY optimized scenarios the observed 
excess was only 8 events and the expected efficiency was slightly lower but 
this scenario required a large transverse momentum imbalance which presents 
difficulties for a symmetric decay. If we consider that the standard Higgs pro- 
duction cross section at Tevatron is over 10^ fb, we see that if the generational 
LLE couplings were the same we would have expected at least a 40 event ex- 
cess in like sign dilepton events; as two neutralinos decaying through a virtual 
charged slepton will produce like sign dileptons half the time. Taus were not 
covered in the dilepton search since they looked for well isolated lepton pairs 
that could be tracked back to a single vertex. This tells us that the LNV 
operator must favor third generation processes in order to avoid constraints. 

Finally, since these decays involve virtual sfermions and small couplings 



5 



they may not be prompt. The general decay length has the form 



L ~ -T^f^^ (3-1) 

where A is the RPV coupling, rrirh is the sfermion mass and N is a factor 
specific to the decay. If the decay length is long enough, the decay will have 
two displaced vertices that are separated from the primary decay vertex. 
However if the decay length is too long, the decay will occur outside of the 
detector and will thus be ruled out by existing missing energy searches. In 
this way we may put bounds on the RPV coupling and sfermion masses. We 
go on to do specific analyses of gaugino decays. 

3.2 Neutralino LSP 

The possibility of a light neutralino LSP decaying via BNV was discussed in 
detail in ref |16j . Here we discuss the possibility of a light neutralino decaying 
via LNV as well. 

The LLE signal would look like Higgs to 4 leptons plus missing energy, 
where the neutralino decays through a virtual lepton. For the LQD channel 
the neutralino may decay either through a virtual slepton or squark. The 
signal would be Higgs to missing energy plus 4 jets, 2 leptons plus four jets, 
or 4 jets one lepton plus missing energy. All of the signals for LNV decay of 
neutralinos may have multiple b squarks. 

BNV decays of the Higgs all have a 6 jet topology, however, the flavor 
antisymmetry of BNV operators does not allow two down type squarks of 
the same flavor to be produced by the same vertex. 

Since the decay happens through a virtual heavy squarks or sleptons, 
we may calculate the decay length. Here we calculate the decay length for 
neutralinos decaying via the BNV operator. 

L - '-^^^:^'4(^'') (3.2) 

a\U2i\ m^x 

5 



3/im nO-^\ / TTifn Y /30 GeVy 



where It/21 1 is an element of the mixing matrix and p-^ is the neutralino's 
momentum. The neutralino may decay via BNV or either LNV operator. For 
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Figure 1: Plot of mu term vs tanb for a Higgs to neutralino cross section 
which beats Higgs to bbar by a factor of 5 and satisfies Z width, b to s 
gamma and chargino mass constraints 

large RPV couplings and small offshell masses this decay is prompt. However 
if (A^/m^o)'* is too small the decay will not happen inside the detector and 
this process will be ruled out by searches for invisible decays of the Higgs. 
There is an intermediate range where the decay happens inside the detector 
but there is a secondary displaced vertex. In this case, the decay products of 
each Xo will track back to these vertices, which are separate from the initial 
vertex where the xo's are made. To illustrate the point we have constructed 
a plot assuming the BNV decay of the neutralinos which shows us the RPV 
coupling vs. offshell squark mass for different displacements of the secondary 
vertex. Though it is easy to be in a region of this parameter space with a 
prompt decay, there is also quite a bit of space that allows for a displaced 
vertex of a few microns. 

Some attention should be paid to the specific issue of like sign dilepton 
constraints in LNV scenarios. As mentioned earlier, like sign dilepton mea- 
surements force us to consider only decays with taus leptons in the final state. 
For decays proceeding through the LQD operator we may easily choose to 
turn on RPV couplings such that only a tau will appear in the final state. 
The LLE operator requires each neutralino decay to two charged leptons and 
a neutrino. Since the LLE operator is antisymmetric in left handed fields, 
and we require the charges leptons to be taus, the neutrino may not be a 
tau neutrino. For example the decay may proceed via the coupling 323 to 
two taus and a muon neutrino. In this case the coupling 233 has an equal 
magnitude and the decay products are just as likely to contain a tau a muon 
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Figure 3: For neutralino masses of 25 and 35 GeV, squark mass vs RPV 
coupling for displaced vertex of 4 microns, 4cm, and Im 

and a tau neutrino. In this case fracl8 of the events will contain like sign 
(non-tau) dileptons, and if we a standard Higgs production cross section and 
an efficiency of 8 percent at 1 fb~^ we expect to see 10 events over back- 
ground at 95 percent confidence. This is within the inclusive excess for the 
Tevatron search. 

We now mention constraints and discuss the parameter space of these 
decays. 

In the case of UDD and LQD operators, there is a contribution of the 
Z width to hadrons and the LLE operator contributes to the total Z width. 
Since our lightest neutralino is mostly Bino, we may sufficiently decouple 
from the Z to suppress large contributions to the width. All scans show 
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Figure 4: Plot of mu term vs tanb for a Higgs to gluino cross section which 
beats Higgs to bbar by a factor of 5 and satisfies Z width, b to s gamma and 
chargino mass constraints 

points within Icr of the measured value. We require that the decay of Higgs 
to neutralinos beat Higgs — > 66 by a factor of 5. In addition we have used 
to decoupling of the Ml parameter from M2 to ensure that we may satisfy 
the chargino lower bound of 102.7 GeV. We have also plotted points such 
that the value of 6 — S7 does not exceed 2a of the measured value while 
scanning over stop mass parameters and charged Higgs masses. We have 
plotted points scanning for a Higgs mass of 91 GeV, however small regions 
of parameter space exist for Higgs masses as low as 87 GeV. 

3.3 Gluino LSP 

The second possibility for gaugino LSPs are gluinos. Light gluinos may eas- 
ily dominate the Higgs decay width as shown in [T7|. Again, gluinos less 
than half of the Higgs mass are maximally incompatible with the standard 
gaugino mass ratio predictions. However light gluinos do have a model build- 
ing benefit, light gluinos do not contribute substantially to squark masses in 
loops. When gluinos have mass, there is a two loop gauge mediation-like 
contribution to squark masses which goes like ^M^2log{^). Here A is the 
cutoff of the theory, and M3 is the gluino mass parameter. If the cutoff of 
the theory is high the squarks get a large mass contribution and the Higgs 
sector may become tuned as noted in [T8] . 

The Higgs decays to gluinos through a triangle diagram involving top 
squarks and top quark. The diagram has a loop factor and is cut off by 
the largest scale in the loop. A similar diagram generates Higgs to gluons 
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except only tops or bottoms appear in the loop. In this case however, the 
coupling of Higgs to tops is only the Yukawa couphng which is not enhanced 
by mixing. The stop's left right eigenstates are not their mass eigenstates. If 
the left right eigenstates have a large mixing, one stop will be much lighter 
than the other. In addition the Higgs coupling to titi is proportional to 
the off-diagonal element of the mass matrix ~ m^^^ — mj2- As the coupling 
becomes large, this decay may easily beat Higgs — > bb. See figure H] for the 
stop sector parameters needed to beat the bbar cross section by a factor of 5, 
where we have plotted parameter space in the limit of light gluinos ~ 20GeV 
or less and a 100 GeV Higgs. 

The gluinos are then free to decay via an LNV or BNV operator. The 
topology of this decay is similar to that of the neutralino, the decay will 
proceed trough a virtual squark, and result in six fermions in the final state. 
In the case of BNV, the gluino will decay to a quark and an offshell squark, 
which will then decay to two squarks; thus the signal is Higgs to six jets. 
If the gluino were to decay via an LNV vertex, it would be via the LQD 
operator. In this case the gluino would decay to a quark and an offshell 
squark. The squark then decays to a quark and lepton. The signal would 
then be Higgs to 4 jets plus missing energy or Higgs to 4 jets plus 2 charged 
leptons. Because there is no flavor antisymmetry in the LNV operator, these 
signals may contain multiple quarks of heavy flavor. In particular, if the 
final state leptons are neutrinos, we may have the signals 4b plus missing 
energy, 4c plus missing energy or 2b and 2 c plus missing energy(since each 
neutralino is free to decay to bottom or charm pairs). Perhaps the most 
striking signal would contain 2b, 2c and 2r. There may also be signals in 
which one gluino decays with a neutrino in the final state, while one decays 
with a tau. In this case the decay products may contain an odd number of a 
heavy quark such as 3b, c, tau, and missing energy or 3c, b, tau and missing 
energy. 

One might worry that a similar process also contributes to Z to gluinos 
which would negatively effect the measured e~^e^ hadrons, and the total 
hadronic width of the Z. However the coupling of stops to the Z is only a 
gauge coupling which can't compete with the coupling to the Higgs. Further, 
the coupling of the stops or sbottoms (but not both at the same time) to 
the Z may be tuned away while remaining large for the Higgs. The b s'j 
measurement mostly constrains allowed values of tanj3, as part of the b ^ s'j 
amplitude goes like l/cos/3. Our plot shows parameter space for which the 
b ^ sj prediction is within 2a of the measured value. Like the neutralino. 
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since the decay of the LSP proceeds through an offshell sparticle, there exists 
the possibihty of a secondary displaced vertex. For sufficiently large RPV 
coupling and small mass of the offshell sparticle, the decay will occur inside 
the detector. In addition, we note that the decay length in this scenario will 
be more likely to be shorter than that of neutralino LSPs, as it goes like -. 
Finally, as is the case with neutralinos, like sign dilepton constraints require 
that the LNV operator must have a structure that favors third generation 
couplings. 

4 Scalar s 

4.1 Topology and General Constraints 

The Higgs may decay to a pair of scalars, which then each decay directly 
through an R parity violating operator. These decays of the Higgs are usu- 
ally prompt, and have a 4 particle final state topology. These scenarios are 
constrained by the chargino lower mass bound and 6 ^ 57 as were the gaug- 
ino LSP scenario. Since the scalar LSP pair must have opposite charges, like 
sign dileptons do not constrain this scenario. Again we must insure that these 
decays can beat Higgs — > 66 by a suitable factor. However, as we explain, 
mass bounds on the scalar LSPs themselves are of the biggest concern. 

Several experiments have put lower mass bounds or production cross sec- 
tion limits on scalars decaying through RPV. The least stringent bounds for 
f decay are from L3 and OPAL respectively and are ruf > 61 for LLE and 
nif > 74 for LQD couplings. For sbottoms the least stringent limits quoted 
by L3 for UDD decays and is > 55. These limits would seem to exclude 
sbottom and stau as LSPs in our scenario. However none of these searches 
looked for RPV decays of these particles in mass windows much below half 
of the Z mass. Table 1 compiles a list of the lowest masses used in direct 
scalar searches [6j,[7j (table 3),y^,^20j. It was assumed that particles less 
than half of the Z mass could be ruled out with Z width measurements and 
with contributions to e+e~ hadrons. 

As we shall elaborate, third generation scalars may have substantial mix- 
ing in their mass matrices. By varying the scalar mixing angle, we may 
change their coupling to the Z. By choosing a small coupling to the Z, we 
may suppress contributions to e'^e" — + hadrons and the Z width for scalars 
with masses below half of the Higgs mass. In this way we may avoid RPV 
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lower mass bounds for third generation scalars. 

For scalars, the mass eigenstates and left right eigenstates are not the 
same. We may write, 

si = SinOsL + CosOsr (4.1) 
s2 = CosOsl - SinOsR (4.2) 

where 9 is the si S2 mixing angle. The mass matrix for third generation 

scalars is given by, 

/ M2^ + ms2 + D-^ rUsA \ 

\ rUsA M^^ + + L>-^ ) 

where A = Asu — /icotp for up- type, A — Asb — iitan(3 for down- type. The 

mixing angle is 

sin2e = 2mbA/{ml^ - ml^) (4.3) 

For third generation sparticles it is possible to get large off-diagonal terms 
and we may have one mass eigenstate which is much lighter than the other. 
In general, it is not too hard to achieve a hght mass eigenstate less than half 
of the Higgs mass. In addition the couplings of the light mass eigenstate to 
the Z are controlled by the mixing angle, 

9ZS2S2 = gihsinl - Q2sin'^9) (4.4) 

which may be adjusted to vanish. Due to differing charges, the angle at 
which the scalar completely decouples will be different for stops, sbottoms 
and staus. It is worth noting that though third generation scalars can be 
decoupled from the Z, nothing will decouple them from photons. Thus while 
we can suppress the contribution of scalars to e~^e~ hadrons, we cannot 
eliminate it. Henceforward we will concentrate on Higgs decaying through 
stau and sbottom LSPs. 

4.2 Sbottom LSP 

The bottom squark is a good candidate for a light LSP. Its decay may proceed 
either through LNV or BNV operator. In the case of BNV the decay proceeds 
as Higgs to hght sbottoms with each light sbottom decaying to an up and 
down quark directly through the RPV coupling. Because of the down type 
flavor antisymmetry no bottoms appear in the final state. The signal is 
thus Higgs to 4 jets, at most 2 of heavy flavor (charm). LNV decays proceed 
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EXP 


LLE 


LQD 


UDD 


DELPHI 


r > 45 




b > 45 


OPAL 


r > 45 


r > 45 


b > 45 


L3 


r > 70 




b > 30 


ALEPH 


r > 45 


T >40, b> 30 


b > 45 



Table 1: LEP2 experiment searches for scalars decaying directly through 
RPV 




Figure 5: Plot of mu term vs tanb for a Higgs to sbottom cross section for 
sbottoms between 7.5 and 30 GeV which beats Higgs to bbar by a factor of 
5 and satisfies Z width, b to s gamma and chargino mass constraints 

through the LQD operator, with the signal two quarks and two leptons. In 
the case of the LQD operator there is no flavor antisymmetry, thus there may 
be a b and a 6 in the final state. If the final state leptons are neutrinos, the 
most visible decay will be Higgs to 2b plus missing energy. The most visible 
signal in this case would be /i — 2c + 2r. 

We will now discuss constraints and the parameter space. Long lived 
sbottoms are ruled out under 92 GeV. However, as long as the RPV couplings 
are large enough and the sbottom is sufficiently more massive than its decay 
products, the decay should be prompt. We see that by adjusting the mixing 
angle of the sbottoms, we decrease the coupling to Z. Following the formula 
quoted in the beginning of this section, we see that this coupling is turned off 
when sinO ~ .39. A lower bound on sbottom masses has been set at 7.5GeV 
by measuring contributions to the over all cross sections of e^e~ — > hadrons 
while turning off the sbottom coupling to Z's [21j. As quoted in table 1, the 
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L3 experiment put an effective upper bound of 30 GeV on sbottom squarks 
decaying through RPV. Previous bounds had been set on the mass sphtting 
of very light sbottoms and the lightest stop |22]. These took into account a 
large stop loop contribution to the Higgs mass and limits on the p parameter 
and required a light stop lighter than 300 GeV. In our scenario however 
the bound is even more relaxed since we do not require a large stop loop 
contribution to the Higgs, and our stop sbottom splitting is not quite as 
extreme. 

The decay rate of Higgs to sbottoms must beat that of bottoms by a 
factor of a few. The ratio of decays rates is given by [23] as 

T-jVb = ptanP^ /2mlsin2e\l - 4m^/m^)^/^ (4.5) 

which gets large for large values of fi and tan/5. In this scenario there are 
5 free parameters; the A term, p, tan/?, and the soft masses. We have 
plotted fi vs tan/? in parameter space for this window of allowed sbottom 
mixing angles. We see in the plot of that lower values of fi and tan/5 are 
ruled out by the upper mass hmit on the sbottom. If the product /xtan/? 
is too small, the off diagonal elements of the mass matrix will not be big 
enough to produce a light mass eigenstate. We therefore expect to find 
parameter space at larger values of tanP than we did for neutralino LSPs. 
The measurements of e~^e~ — *■ hadrons constrain allowed values of 6. Janot 
tells us the largest allowed mass range for sbottoms occurs between mixing 
angles .3 < sinO < .45 [21] (see fig 8). For A terms and fi terms of a few 
hundred GeV we see that we may fall into this range of mixing parameter 
if A and fi cancel to within 20 or 30 percent of their value. The constraint 
from b to s gamma follows as before and we scan over stop sector parameters 
and obey the chargino constraint. However neither of these constraints is as 
restrictive as for the case of light neutralinos since we are free to make the 
gaugino sector heavier in general. 

4.3 Stau LSP 

The decay of the Higgs may proceed through a stau anti-stau LSP pair. The 
direct RPV decays of the tau happen through LNV operators only. 

For LQD decays of the stau, the final state will have 4 quarks, 2 up type 
and 2 down type. Because there is no flavor antisymmetry, the final state 
may consist of all quarks of heavy flavor, h bbcc. 
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Figure 6: Plot of stau mixing vs hadronic cross section at the Z pole. The 
lower region is ruled out by PEP PETRA and TRISTAN measurements, the 
regions to the left and right by Z pole data from LEP 1. 
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Figure 7: Plot of allowed stau mass vs. stau mixing angle 



For decays with quarks in the final state, the f masses below half the Z 
mass are constrained by the measurement of e^e~ —>■ hadrons. Following the 
method of Janot for light sbottoms, we may go back to LEP2, PEP, PETRA, 
and TRISTAN measurements of the total e^e~ hadrons cross sections to 
set lower mass bounds on the staus [21j . Ref [2T] contains a compilation of the 
total ahad measurements for various experimental energies. We first calculate 
the absolute minimum allowed f mass by assuming complete decoupling from 
the Z and calculating the resulting contribution to e~^e~ — >■ hadrons for 
production from photons alone for different experimental energies. We then 
perform a fit for the hadronic cross section contribution. We find that 
at 95 percent confidence the lower limit on stau masses is 11 GeV with a 
best fit value of 57 GeV. A plot of the LEP and low energy cross section 
measurements and cross section predictions for light staus appears in the 
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appendix. There seems to be significant parameter space for light stuas. We 
must note that Janot performed a more sophisticated analysis of the low 
energy data which may further constrain the lighter stau masses. However 
for stau mass heavier than 28 GeV the lowest energy data does not constrain 
our scenario. We will therefore conservatively consider the viable allowed 
stau mass window to be above 28 GeV. 

In addition we may analyze the Z pole data from LEP 1. As quoted 
by Janot, the Z pole measurements limit new contributions to the hadronic 
cross section to 56pb at 95 percent confidence. Taking these measurements 
into account, we may exclude stau masses for different couphngs of staus 
to the Z. When the f has small coupling to the Z, we find a mass window 
between 28 and 45 GeV. Following the formula in section 4, we see that the 
Z coupling to f is turned off when sinO = .67, and allowed mixing angles 
range from .6 < sin6 < .7 in the case of light staus. We have plotted the 
total contribution from light stau production at the Z pole vs mixing angle, 
as well as allowed stau mixing angles vs stau mass taking into account both 
Z pole and PEP,PETRA and TRISTAN data. 
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Figure 8: Plot of mu term vs tanb for a Higgs to staus cross section, with 
a Higgs mass of 105 GeV, which beats Higgs to bbar by a factor of 5 and 
satisfies Z width, b to s gamma and chargino mass, and hadronic cross section 
constraints 

In principle the parameter space for this decay looks much like the pa- 
rameter space for sbottoms. However, the constraint on the allowed mixing 
angle is tighter since a large contribution of e+e~ —>■ hadron comes directly 
from the photon coupling to tau with no coupling to the Z at all. Again get- 
ting a light stau requires large contributions to the off-diagonal mass matrix 
elements, in this case given by mT-{A — fitanjS). Since rrir is small, a large 
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product /i tan(3 is required for large mixing, even larger than what was re- 
quired for light sbottoms. We keep in mind that we do not want excessively 
high values of fi, since this would again introduce tuning into the Higgs po- 
tential. In addition tan/3 must not become too large or the bottom Yukawa 
coupling will become nonperturbative. We also need to choose fi and tan(3 
such that we may obey the lower mass bound on the lightest chargino. In 
this case, the variables effecting the stau sector that are at play in the b —>■ s'j 
process are only tanP and fi, the stau and stop sectors are fairly uncoupled 
and 6 — >■ S7 is not a strong constraint in the light stau scenario. We have 
plotted parameter space for all constraints. 

The LLE decay of the stau is an interesting subject. If allowed, each stau 
would decay into a charged lepton plus a neutrino, and the signal would be 
opposite sign dileptons plus missing energy. The most striking signal here 
would he h ^ 2r+ /E. A search with similar signature, gauge mediated 
decays of stau to tau plus gravitino [23], highly constrains this scenario. 
This search rules out the tau plus missing energy signal for stau mass larger 
than 2GeV. It should be noted that this search found a small window for 
mr < ruf < 2GeV. One might guess that in the case where both staus 
decay hadronically, contributions to e~^e~ hadrons rule out this possibility. 
However, we see that the PEP, PETRA, and TRISTAN data placed cuts 
specifically to rule exclude the background from r"*"r~ and it is likely these 
decays may have been missed. See for example [2S])- Thus there may be a 
signal for Higgs to 2 stau plus missing energy in the improbable event that 
there exist super-light staus. In this case, since there is a only a small window 
between the tau and the stau, if the RPV coupling is small, the stau may live 
for some time before it decays. In this case, in addition to a 2r plus missing 
energy decay, they may also be a secondary displaced vertex. For example 
for LLE coupling of the size a few times 10~^ a 1.8 GeV stau may live for 
100 microns. Stau pair decays with one or two light leptons may also occur. 

5 Conclusions 

In the R parity violating MSSM without standard gaugino mass relationships, 
we find a menagerie of new decays for a light Higgs. We present the list of 
possible LSPs and topologies in the table found above. 

One may wonder how all of these signals pass current bounds. There 
are several current searches that one must consider and we will first mention 
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general of these. DELPHI did an analysis for Higgs goes to anything which 
could put a lower mass bound on the Higgs of 82 GeV There has been 
a Higgs to missing energy search which put a lower bound on the Higgs mass 
of 114 GeV [27]. This is a concern for us only when the Higgs decay is 
not prompt as is the case for gaugino LSPs. It is easy to avoid falling into 
this search if RPV couplings are large or if squark or slepton masses are of 
reasonable size. In fact, in the case of squarks, we would not want masses 
to go much beyond ITeV, for then we encounter the same tuning problem in 
the Higgs potential that we wished to avoid. 

We shall now consider searches which apply to decays with all hadronic 
final states. First there is the 2 jet flavorless search which puts a lower 
bound on the Higgs of 113 GeV [2B]. In order for this search to be sensitive 
to our scenario, we would have to force our 4 and 6 jet final states into two 
jets. Such forcing is used in a set of searches where members of the Higgs 
multiplet undergoes cascade decay [29] [30]. Thus we worry about decays 
e^e~ — > H2Z HiHiZ, where Hi decays to bottoms or taus. In the all 
hadronic case, the final state is 4b. The event is forced into two 2b jets 
and the efficiency is calculated of a 4b even being picked up by a 2b search. 
Using the DELPHI search table we can guess our efficiency of being picked 
up by the 2b search. This should give us a good idea about matching our 
4 jet signals to 2 jets, for 6 jets the efficiency would be even worse. As an 
example, notice that to rule out the Higgs at 80 GeV, the 4b search needs 
about 5.5 times the number of events as the 2b. Even assuming both of our 
b's were tagged, our scenario would only be picked up 18 percent of the time. 
The flavorless search is the same as the 2b search without the b tags|28]. 
LEP performed a search in which HZ was produced and Higgs decayed to 
WW* [32] • An analysis was done for the final state Z —>■ uu, H —>■ qqqq, 
which is relevant to our light stau scenario. This search places a lower mass 
bound on the Higgs of 105 GeV. In addition the WW* search considered a 
final state in which e'^e~ HZ — > Qq. In this case however, cuts are applied 
which reconstruct the masses of the final state Z and Ws, so this search is 
not likely to be sensitive to our 4 or 6 jet signals. 

We now turn to searches with quarks and charged leptons in the final 
states. The cascade decays mentioned above fall into this category produce 
final states with 4r and 2b+2r. We can compare these final states to all 
of ours and see that these searches are not directly sensitive to any of our 
final states. Of our final states without missing energy these signals come 
closest are the six body decay 2b+2c+2r and the four body decay, 2c+2r. 
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Reconstruction of the cascade decay requires multiple b tags. It is unlikely 
that the 2 charms will receive a large b-likeness parameter. In addition to the 
b tag problem, our 2b+2c+2r decay must be force four jets into two loosing 
efficiency. We do not expect that these searches will constrain our scenario. 
Cascade decays with even more final state b's and rs occur through the 
process e"'"e~ H2H2 HiHiHi. Again, none of these final states match 
directly to our signals. Events with 4 or 6 b's require even more b tags, which 
are unlikely to match up to our 6 particle final states that have at most 2 b's 
and none of our final states contain more than two taus. 

Finally there are decays with missing energy in the final state. A com- 
pendium of some searches that constrain Higgs signals with missing energy 
is found in [3T]. In particular this paper quotes upper bounds on the Higgs 
production cross section for Higgs decays with final states and 2q+ /E. These 
limits do not come from a Higgs search but from the LEP2 squark searches 
where sparticles are pair produced and decay to quarks and neutralinos. This 
signal is identical to our scenario with quarks and missing energy in the fi- 
nal state. In this lower bound may be set on the Higgs mass of 103 
GeV if the final state quarks are light and 111 GeV if both of the final state 
quarks are b quarks. These bounds do not apply in cases where the final 
state quarks are not the same fiavor. One might consider the cascade decays 
e^e~ — >■ -^2-2^ ^ HiHiZ Ab/Ar Z. If in this case the Z decays invisibly the 
overall signal would be identical to our H — > 46/4r + 2z/ final state. The 4b 
search involved cuts which considered parameters such as the missing mass 
and the ln{xmz) parameter where the missing mass is forced to the Z mass. 
These cuts are bound to exclude data in our scenario. Further we can see from 
the exclusion plot in fig 12 that if the efficiency of observing out scenario is 
60 percent, all of parameter space is open up down to the 82 GeV Higgs mass 
bound. The 4r search constraints on the decay topology such as restrictions 
on the angle between charged particle pairs, in addition to cuts on the miss- 
ing momentum which make it insensitive to our decay. However the in the 
WW* search an analysis was done for the channel e^e~ HZ WW*Z 
where the Ws decay hadronically and the Z invisibly. Here we cannot avoid 
constraints and a lower limit of 105 GeV may be placed on the Higgs mass. 
Finally, the WW* search analyzed a Higgs channel with charged light leptons 
and missing energy in the final state. This signal is identical to out final state 
where scalar LSPs decay through the LLE operator. Even in the case that 
we consider an LLE operator like 313, where the final state leptons may be 
only taus, the antisymmetry of the RPV operator forces us the coupling 133 
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LSP 


LLE 


LQD 


UDD 


Xo 


4r+2z/ 


4b/4c+2i/, 2b+2c+2i/, 2b+2c+2r,3b+c+r+z/, b+3c+r+z/ 


2b+2c+2q 


g 




4b/4c+2i/, 2b+2c+2z/, 2b+2c+2r,3b+c+r+z/, b+3c+r+z/ 


2b+2c+2q 


b 




2b+2i/, 2c+2r, b+c+z/+r 


2c+2q 


r 


2r+2z/ 


2b+2c 





Table 2: Higgs decay signals for all possible LSPs and RPV operators 



to be the same size. Thus we get light leptons a quarter of the time and this 
search will also constrain the 2 tau plus missing energy channel. The search 
places an upper bound on the Higgs production cross section of .044pb and in 
our case, this translates to a Higgs mass lower bound of around 104 GeVfor 
final states with taus or light leptons plus missing energy. We may lower the 
bound by turning on more than one RPV coupling at once and decreasing 
the likelihood that the decay products are symmetric. For example, turning 
on two coupling at once which are equal in size decreases the Higgs mass 
bound to 95 GeV. 

Other asymmetric decays are also listed in table 2 which involve 2 RPV 
coupling being turned on at once. These scenarios are bizarre enough to be 
mostly unconstrained. The final state b+c+r + z/ is mentioned in the WW* 
search, however the combinatorics of turning on two RPV operators at once 
allow us to avoid the Higgs mass bound of 95 GeV set by this scenario. 
Stranger 6 body decays listed in our table have are not constrained by any 
relevant search. 

Overall we have proposed over a dozen distinct channels for Higgs discov- 
ery with multiple heavy flavor particles in their final states. Some of these 
decays would be very hard to detect, for example those decays in which the 
final state is 4 or 6 jets. However, many of these decays have interesting 
missing energy signatures, some of which are quite bizarre - for example the 
b+3c+r +u decay of the gauginos. Those decays that proceed through a 
gaugino LSP have the added bonus of possible secondary displaced vertices. 
We may imagine modifying existing searches to look in some of these chan- 
nels, for example modifying existing 2 or 4 b and r searches to be sensitive 
to missing energy. In addition we might hope to detect events which contain 
bottoms at LHCb as has been recently proposed [33] . 
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LSP 


Signature 


Mass Bound 


Search 


Xo 


Aq + 2v 


105 GeV 


WW* with invisible Z decay 


9 


Aq + 2v 


105 GeV 


WW* with invisible Z decay 


b 


2q + 2u 


103 GeV 


SUSY squark search 




2h + 2v 


111 GeV 


SUSY squark search 




Aq 


105 GeV 


WW* with invisible Z decay 


f 


TT + 2U 


104 GeV 


WW* 




ri + 2v 


104 GeV 


WW* 




4q 


105 GeV 


WW* with invisible Z decay 



Table 3: Higgs Mass Lower Bounds for Various Channels. For decays not 
listed current searches do not severely constrain the Higgs mass. 



6 Appendix 

Below we have plotted hadronic cross section vs center of mass energies. The 
points are the measurement of the hadronic crossection minus the SM theory 
prediction with two sigma error bars. The curves are contributions to the 
hadronic cross section from the decay of low energy staus of various masses. 
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Figure 9: Plot of Energy vs hadronic cross section for low energy data and 
stau contributions, from the top down the curves are for staues of mass 10, 
15, 22, 28, and 45 GeV 
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